function [ error ] = compareasymptotic(N, sigma)

D = 1;

epsilon = 2*pi*sigma/N;

[r, theta, v] = equalspacepolar(N,sigma);

va = zeros(length(r),length(theta));
x = va;
y = va;

va0 = -r.^2./(4*D) + 1/(4*D) - 1/(2*D*pi*sigma) * epsilon * log(sin(pi*sigma/2));

for i = 1:length(theta)
   va(:,i) = va0; 
end

for i = 1:length(r)
    for j = 1:length(theta)
        x(i,j) = r(i)*cos(theta(j));
        y(i,j) = r(i)*sin(theta(j));
    end
end

error = abs(v-va);
%surf(x,y,error,'Edgecolor','none');
contourf(x,y,error, [min(error(:)):(max(error(:))-min(error(:)))/30:max(error(:))],'Edgecolor','none');
colorbar;
end